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PROGRAM  SPECIFICATIONS: 

Language :  ANSI  FORTRAN  (FORTRAN  IV) 

Solution  Requirements :  The  run  command 

RUN  WESLIB/CORPS/H612U,  R 

plus  the  input  variables  defined  in  (e) 

Method  of  Analysis :  Solves  algebraic  equations  for  the  specific 
energy  and  specific  force. 

Core  Requirements  G635:  10  K  words 

External  Storage:  None 

Restrictions :  Velocity  distribution  is  uniform  across  the  flow 
section  and  the  invert  slope  is  small  (i.e. ,  <10  degrees). 

General  Equations: 

Q2 

Specific  Energy  E  (ft  of  water)  E  =  y  +  - — 

2gA 

q  _  02 

Specific  Force  F  (ft  of  water)  F  =  Ay  +  — — 

where:  Q  is  the  discharge  (cfs),  g  is  acceleration  of  gravity 
2 

(32.2  ft/sec  ),  y  is  flow  depth  (ft),  A  is  cross-section  area 
(ft  ),  and  y  is  the  distance  (ft)  of  the  centroid  of  the  area  A 
below  the  water  surface. 

Range  of  Quantities:  Unlimited  for  practical  application. 

Accuracy:  Governed  by  accuracy  of  input  data;  the  specific  energy 
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and  specific  force  are  printed  to  +0.01  ft  and  ft  . 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  Computer 
Program  Library  Standards  and  Documentation,  Appendix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  722-F3-R0-6AJ 

2.  TITLE:  H6124  -  Specific  Energy  and  Specific  Force  in  Either  a 
Trapezoidal,  Triangular,  or  Rectangular  Open  Channel. 

3.  REVISION  LOG:  N/A 

4.  PURPOSE  OF  PROGRAM:  To  compute  the  specific  energy  and  specific 
force  in  either  a  trapezoidal,  triangular,  or  rectangular  open  channel. 
Reference : 

Chow,  Ven  Te,  Open  Channel  Hydraulics,  McGraw-Hill,  1959,  PP  54. 

5.  STEP  SOLUTION: 


a.  Enter  inputs: 


(1)  Discharge  (Q),  cfs 

(2)  Channel  bottom  width  (b),  ft  (zero  if  triangular  section) 

(3)  Channel  side  slope  (Z)  expressed  as  the  cotangent  (zero  if 
rectangular  section) 

(4)  Flow  depth  (y),  ft 

b.  The  computational  steps  are: 

,  2 


(1)  Area  (A),  ftc 


A  =  (b  +  yZ)y 


(2)  Specific  Energy  (E),  ft  of  water 

Q2 

e  =  y  + 

2gA^ 


where  g  is  acceleration  of  gravity  (32.2  ft/sec  ). 
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—  3 

(3)  Moment  (Ay),  ft  ,  of  area  with  respect  to  the  surface 

2  3 

Ay  =  aLb  +  zfz 

W  2  3 

(4)  Specific  force  (F),  ft3  of  water 

F  =  Ay  + 

gA 

c.  The  given  data  and  specific  energy  and  specific  force  are 
printed. 

6.  ACCURACY :  Governed  by  accuracy  of  input  data.  The  specific  energy 
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and  the  specific  force  are  printed  to  jjD.Ol  ft  and  ft  . 

7.  REMARKS:  Velocity  distribution  is  uniform  across  flow  section  and 
the  invert  slope  is  small  (i.e.,  <10  degrees). 


PART  II:  COMPUTER  FUNCTIONAL  DESCRIPTION 
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3.  EQUIPMENT  AND  OPERATING  SYSTEM:  The  program  was  developed  on  a 
G635  time-share  system  in  which  input/output  equipment  consisted  of  a 
Model  33  remote  teletype.  It  is  now  operational  on  the  WES  G635, 
Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 

4.  INPUT  REQUIREMENTS:  The  required  inputs  are  entered  via  the  user's 
time-share  terminal  device  in  free  format.  All  input  cues  and  reads  are 
performed  in  the  main  program.  The  subroutine  handles  the  computations. 
Since  computations  are  done  in  the  subroutine,  the  necessary  inputs  to 
subroutine  H6l24  are  passed  via  the  CALL  statement.  The  calling 
sequence  is: 

CALL  H6l2U(arg1, . . . ,argg) 

where : 

arg^  -  discharge,  cfs 

arg2  -  bottom  width,  ft 

arg^  -  side  slope,  cotangent 

arg^  -  flow  depth,  ft 

arg<_  -  specific  energy,  ft  of  water 

3 

argg  -  specific  force,  ft  of  water 
All  values  are  floating  point.  Arguments  1-4  are  inputs  and  arguments 
5  and  6  are  output. 

5 .  SECONDARY  STORAGE  INPUT:  None 

6.  INPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the  input 
variables  in  program  H6124. 
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DISCII  -  discharge,  cfs  (srq) 
BWID  -  channel  bottom  width,  ft  (argp) 
SSLOPE  -  channel  side  slope,  cotangent  (arg.J 
NDEPTHS  -  number  of  flow  depths,  integer 

Y  -  flow  depth,  ft  (arg^) 


7.  OUTPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the  out¬ 
put  variables  in  program  H6l24: 

ENERGY  -  specific  energy,  ft  of  water  (arg. ) 

FORCE  -  specific  force,  ft^  of  water  (arg(0 

8 .  PROGRAM  ERROR  MESSAGES :  None 

9.  VARIABLE  DEFINITIONS: 

2 

AREA  -  flow  depth  area,  ft 

3 

AY  -  moment  of  area  with  respect  to  surface,  ft 

BWID  -  channel  bottom  width,  ft 

NDEPTHS  -  number  of  flow  depths,  0  <  NDEPTHS  £  25 
Y  -  flow  depth,  ft;  dimensioned  max  25 

DEPTH  -  flow  depth,  ft;  equal  to  Y(l)  for  I  =  1 , . . . , NDEPTHS 

DISCII  -  discharge,  cfs 

ENERGY  -  specific  energy,  ft  of  water 

FORCE  -  specific  force,  ft^  of  water 

2 

G  -  acceleration  of  gravity,  32.2  ft/sec 

HFILE  -  five  character  name  of  program;  passed  to  WESLIB  count- 
routine  HACCT 

JKL  -  direct  return  from  RERUN  to  desired  input  read 
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KKK  -  total  number  of  inputs,  passed  to  RERUN 

LQX  -  equal  1,  print  instructions  from  RERUN ;  equal  3,  no 

print 

LQZ  -  equal  1,  execute  all  input  cues  and  reads ;  equal  2,  call 

WESLIB  routine  RERUN  and  enter  only  desired  inputs 

SSLOPE  -  channel  side  slope,  cotangent 

TYPE  -  character;  type  of  cross  section,  either  trap,  rect,  or 
tria. 

ZZZZZ  -  character;  equal  RE,  rerun;  equal  ST,  stop 
10.  EXAMPLE  CASE:  Compute  the  specific  energy  and  specific  force  for 
6  given  flow  depths  of  a  rectangular  channel  cross  section, 
a.  Input  data: 

Bottom  width  (0  if  tria  sect)  (BWID)  =  3.51  ft 
Side  slope  (0  if  rect  sect)  (SSLOPE)  =  0 

Discharge  (DISCH)  -  100  cfs 
Number  of  depths  (NDEPTH)  =  6 

Flow  depth(s)  (Y(I)  for  1  =  1,  NDEPTHS )  =  1,  1.1*3,  2.5,  U,  5.27,  6.5  ft 
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b.  Output : 


INPUT  H6124  -  SPECIFIC  ENERGY  AND  FORCE-TRAP , RECT , TRIA-OPEN  CHANNEL 

AA-ENTER  CHANNEL  BOTTOM  WIDTH  IN  FT (ZERO  IF  TRIA  SECT). 

=  3.51 

AB-ENTER  CHANNEL  SIDE  SLOPEIZERO  IF  RECT  SECT). 

=  0 

AC-ENTER  DISCHARGE  IN  CFS. 

=  100 

AD-ENTER  THE  NUMBER  OF  DEPTH! S )  FOR  WHICH  SPECIFIC  FORCE  AND  ENERGY  ARE 
TO  BE  CALCULATED.  MUST  NOT  EXCEED  25  DEPTHS. 

=  6 

AE-ENTER  THE  6  DEPTH! 5 )  SEPARATED  BY  COMMAS. 

=1,1.43,2.5,4,5.27,6.5 


OUTPUT  H6124  -  SPECIFIC  ENERGY  AND  FORCE  OF  A  RECT.  OPEN  CHANNEL 

BOTTOM  WIDTH  =  3.51  FT. 

SIDE  SLOPE  =  0.  H :  1 V 

DISCHARGE  =  100.00  CFS. 


DEPTH 

SPECIFIC 

SPECIFIC 

ENERGY 

FORCE 

!  FT ) 

!  FT  ) 

( FT  *  *  3 ) 

1.00 

13.60 

90.23 

1.43 

7.59 

65.46 

2.50 

4.52 

46.36 

4.00 

4.79 

50.20 

5.27 

5.72 

65.53 

6.50 

6.80 

87.76 

ENTER  RERUN  OR  STOP 
=  STOP 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  III:  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE:  H6124  -  Specific  Energy  and  Specific  Force  in  Either  a 
Trapezoidal,  Triangular,  or  Rectangular  Open  Channel 

3.  SOURCE  PROGRAM  LISTING:  See  pages  9-11. 

U.  NUMERICAL  AND  LOGICAL  ANALYSIS:  No  special  numerical  techniques 
used. 

5 .  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  None 

6.  MISCELLANEOUS :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational  on  the  WES 
G635,  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 
The  source  listing  on  page  9  contains  the  first  line  run  command  and 
brief  for  E6l2k.  This  first  line  run  command  runs  the  binary  H612LB 
of  the  source  listing  on  pages  10-11  (Fortran  source  of  H6l2^)  and 


attaches  the  WESLIB  routines  HACCT  and  RERUN. 


H  612^4 


000 1  * ffRUN  WESLI8/C0RPS/H6124B,R;WESLIB/RERUN,R;WESLIB/HACCT ,R 
0800  63THIS  PROGRAM  COMPUTES  THE  SPECIFIC  FORCE  AND  ENERGY  FOR  A  TRAP- 
0810  62  EZO I  DAL  »  TRIANGULAR,  OR  RECTANGULAR  CROSS-SECTION.  INPUTS  RE- 
0820  6 1QUIRED  ARE  BOTTOM  WIDTH  IN  FT (ZERO  IF  TRIANGULAR),  SIDE  SLOPE 
0830  6 1  EXPRESSED  AS  THE  COTANGENT  OF  THE  ACUTE  ANGLE  WITH  RESPECT  TO 
0040  63THE  HOR I ZON T A L ( ZERO  IF  RECTANGULAR ) ,  DISCHARGE  IN  CFS,  THE  NUM- 
0850  63BER  OF  DEPTHS  FOR  WHICH  THE  SPECIFIC  FORCE  AND  ENERGY  ARE  TO  BE 
0360  60COMPUTED,  AND  THE  DEPTH(S).  OUTPUT  INCLUDES  THE  GIVEN  DATA, 
0870  62TYPE  OF  CROSS-SECTION,  AND  A  TABLE  OF  DEPTH! S )  AND  CORRESPOND- 
0880  30ING  SPECIFIC  FORCE  AND  ENERGY. 

0999X06FINISH 
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00001*#RUN  *=;/C0RPS/H6124B(N0G0) 

10000  CHARACTER  T YPE*4 , HF I L E#5 
10010  DIMENSION  Y(25) 

10020  HFI L  E ;5MH6 124 

10030  LQX-  1 ; LQZ-  1 

10040  15000  CALL  HACCT (HFILE) 

10050  PRINT  15111 

10060  15111  FORMAT ( /"  I NPUT  H6124  -  SPECIFIC  ENERGY  AND  FORCE-TRAP , RECT , 
100708TRIA-OPEN  CHANNEL"//) 

10080  GO  T0(10,15),LQZ 
10090  10  PRINT  20 

10100  20  FORMAT (48HAA-ENTER  CHANNEL  BOTTOM  WIDTH  IN  FT ( ZERO  IF  TRIA 
1011087H  SECT) . ) 

10120  25  READ , BWID 
10130  GO  TO(30,15),LQZ 
10140  30  PRINT  35 

10150  35  FORMAT (47HAB-ENTER  CHANNEL  SIDE  SLOPECZERO  IF  RECT  SECT).) 

1  0160  40  READ , SSL  OPE 
10170  GO  TOC45, 15) , LOZ 
10180  45  PRINT  50 

10190  50  FO RM AT (26HAC- ENTER  DISCHARGE  IN  CFS.) 

10200  55  READ ■ DISCH 
10210  GO  TO(60,15),LQZ 
10220  60  PRINT  65 

10230  65  FORMAT ( 4 9HAD-ENTER  THE  NUMBER  OF  DEPTHCS)  FOR  WHICH  SPECIFI 
10240822HC  FORCE  AND  ENERGY  ARE/33HT 0  BE  CALCULATED.  MUST  NOT  EXCEE 
1025081 2HD  25  DEPTHS. ) 

10260  70  READ . NDEP  THS 

10270  I F ( NDEPT  HS-25  )  75,75.80 

10280  80  PRINT ."NUMBER  DEPTHS  >  25, RE-ENTER  NUMBER  <  OR  =  25" 

10290  GO  TO  70 

10300  75  GO  TO(85,15),L9Z 

10310  85  PRINT  90 , NDEPTHS 

10320  90  FORMATC'AE-ENTER  THE  ",I2,"  DEPTHCS)  SEPARATED  BY  COMMAS.") 
10330  95  READ. (YCI), 1=1, NDEPTHS) 

10340  GO  TO(100,15),LQZ 
10350  15  KKK=5 

10360  CALL  RERUNCKKK.LQX, JKL ) 

10370  GO  TO(25.40,55,70,95.100),JKL 

10380  100  TYPEr4HTRAP 

10390  IFCBWID)  105.110,105 

10400  110  TYPE-4HTRIA 

10410  105  IFCSSLOPE)  115,120,115 

10420  120  TYPE-4HRECT 

10430  115  PRINT  1 25 , TYPE, BWI D , SSLOPE , DISCH 

10440  125  FORMA TC/46HOUTPUT  H6124  -  SPECIFIC  ENERGY  AND  FORCE  OF  A  , 
104508A4 ,".  OPEN  CHANNEL "//15HB0TT0M  WIDTH  =  ,F9.2,4H  FT. 

1 0 46 0 1/1 3HSI DE  SLOPE  =  , F 1 1 .2,4HH: 1V/12H DISCHARGE  =  ,F12.2,5H  CFS.// 

1 047 08 4X , 5HDEPTH , 2( 4X , 8HSPECI FIC )/14X, 6HENERGY , 6X , 5HFORCE/5X , 4H<  FT ) , 
1048016X.4HCFT) , 6X , 7H ( FTXX3 ) ) 
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10490  DO  140  1=1, NDEPTHS 
10500  DCPIH-Y(I) 

10510  CALL  H6124<0ISCH,BU>ID,SSLCPE, DEPTH, ENERGY. FORCE) 

10520  140  PRINT  1  JO , DEPTH , ENERGY . FORCE 
10550  150  FORMAT ((F9.2,F11.2,F12.2>) 

10540  PRINT,"  " 

10550  102=2 

10560  CHARACTER  ZZZZZ1*Z 

10570  16000  PRINT,  "ENTER  RERUN  OR  STOP" 

10580  READ  16001,  ZZZ7ZZ 
10590  16001  FORMAT  <  A2 ) 

10600  IFIZZZZZZ. EQ.2HRE)  GO  TO  15000 
10610  IFIZZZZZZ. EQ .2HST 1  GO  TO  20000 
10620  PRINT , "ERROR  *#*  RETYPE" 

10650  GO  TO  16000 
10640  20000  STOP ; END 

20000  SUBROUTINE  H6 1 24 < DI SCH . BWI D, SSLOPE , DEPTH , ENERGY , FORCE) 

2C010  G=32.2 

20020  AREA=DEPTH#(BWID+SS10PE#DEPTK) 

20030  ENERGY=DEPTH*DI5CH**2/(2. *G*AREA»*2) 

20  04  0  F0RCE  =  BUIDKDEPTNA«‘2/2  .  +SSLOPEH  DEP  TH*  *  3/3  .  +  DISCH*#2/ ( G#  ARE  A  ) 
20050  RETURN 
20060  END 
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